Spontaneous Ratchet Effect in a Granular Gas 
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The spontaneous clustering of a vibrofluidized granular gas is employed to generate directed 
transport in two different compartmentalized systems: a "granular fountain" in which the transport 
takes the form of convection rolls, and a "granular ratchet" with a spontaneous particle current 
perpendicular to the direction of energy input. In both instances, transport is not due to any 
system-intrinsic anisotropy, but arises as a spontaneous collective symmetry breaking effect of many 
interacting granular particles. The experimental and numerical results are quantitatively accounted 
for within a flux model. 
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The conversion of random fluctuations into directed 
motion in a periodic system - termed ratchet effect - is 
of interest in a wide variety of physical, biological, and 
technological contexts Q. Overall, two necessary condi- 
tions for the emergence of the effect have been identified 
Q , namely (i) the absence of thermal equilibrium and (ii) 
the breaking of the inversion symmetry. The most com- 
mon realization of the latter condition is by means of an 
intrinsic anisotropy in the system, e.g. a spatially peri- 
odic but asymmetric potential, or a so-called dynamical 
asymmetry due to a time-dependent driving force with 
zero mean but non- vanishing higher moments [2( . In this 
Letter, we are concerned with perfectly symmetric sys- 
tems where the symmetry breaking arises spontaneously 
as a collective effect of many interacting particles via an 
ergodicity breaking non-equilibrium phase transition . 
While a proof of principle for this so-called spontaneous 
ratchet effect has been provided theoretically with several 
different "minimal models" |2j,|3|, we explore here for the 
first time a real system, namely a vibrofluidized granular 
gas, experimentally, numerically, and analytically. 

Transport in granular systems due to a ratchet effect 
has been previously demonstrated both for asymmetric 
potentials 4] and for dynamical asymmetries |5|. While 
in both these cases even a single granular particle moves 
in a preferential direction, the spontaneous ratchet ef- 
fect we consider here is a genuine collective phenomenon. 
The underlying mechanism is the enhanced energy loss 
in already dense regions of a granular gas due to the in- 
elastic collisions, resulting in a spontaneous separation 
into dense and dilute regions [6j. This clustering effect 
can be controlled by confining the particles to a series 
of compartments that are connected by small apertures 
(slits) and subjected to gravity |]l|8j. When the granular 
system is fluidized by shaking it vertically, the particles 
cluster into the compartments which are initially or due 
to random fluctuations slightly denser, until a dynamical 
equilibrium is established in which the particle flux out 
of a compartment is balanced by the influx coming from 
its neighbors. 

Hence, all ingredients for a spontaneous ratchet effect 
are present: First, random fluctuations in the form of de- 
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FIG. 1: (a) Schematic cross section of a granular fountain ex- 
hibiting spontaneous symmetry breaking into a "cold" and a 
"hot" compartment and a concomitant spontaneous circular 
particle flow, (b) By folding out the geometry of several ad- 
jacent fountains, and adding cyclical boundary conditions, a 
granular ratchet is obtained with almost unchanged popula- 
tions of the respective compartments and the particle currents 
between them (here shown for K = 4 compartments). 



terministic noise are created intrinsically by the chaotic 
dynamics. Second, the system is kept away from thermal 
equilibrium by the vertical vibrations. Third, in a per- 
fectly symmetric, periodic compartmentalization, stable 
steady states with a spontaneously broken symmetry yet 
still respecting the periodicity seem possible, but are not 
at all obvious: their actual existence is the main finding 
of our present work. 

We start with a granular gas in a container that is di- 
vided into K = 2 equal compartments by a wall with 
a slit at a certain height h. At vigorous shaking, this 
system has a unique steady state with an equal pop- 
ulation of both compartments Q. Upon reducing the 
driving, this symmetric state loses stability via a contin- 
uous phase transition, resulting in a pair of stable so- 
lutions with spontaneously broken symmetry: They are 
clustered states with one diluted and one crowded com- 
partment ■ If we now add a sufficiently small hole at 
the bottom of the wall separating the compartments, the 
clustered states are expected to subsist (see Fig.[TJi), re- 
sulting in a particle flow through the hole, balanced by a 
flow through the slit in the opposite direction. We thus 
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expect a "granular fountain ', i.e. a spontaneous con- 
vective flow, imposed upon the system by the clustering 
phenomenon. 

To verify these predictions, we conducted both experi- 
ments and molecular dynamics simulations. In the exper- 
iments we used N = 407 beads of stainless steel (radius 
r = 1.18 mm, normal restitution coefficient v w 0.9) in 
a high perspex container with a ground surface of 41.5 
x 25.0 mm 2 . When the particles are at rest, this corre- 
sponds to a filling level of 2.0 layers. The box was divided 
into K = 2 equal compartments by a wall, with a horizon- 
tal slit of 25.0 mm by 5.0 mm, starting at height h = 25.0 
mm, and a hole of dimensions 3.0 mm x 3.0 mm at the 
bottom. The entire device was subjected to vertical, sinu- 
soidal vibrations with variable frequency / (ranging from 
47 Hz to 150 Hz) and fixed amplitude a = 1.0 mm. The 
results are given in Fig. [5] in terms of the dimensionless 
control parameter 
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which is motivated ,7] from the kinetic theory for di- 
lute granular gases (see also Eq. J5J below). Here 
(= 19.7 x 25.0 mm 2 ) is the surface area of each compart- 
ment, and g — 9.81 m/s 2 . As anticipated, vigorous shak- 
ing (small B) leads to a symmetric steady state, giving 
way to a clustered state with the predicted fountain effect 
for intermediate driving. For very weak driving (large B) 
the system returns to equipartition via the hole at the 
bottom of the container, since the flux through the slit 
(cf. Fig. QJi) becomes negligibly small. Furthermore, a 
pronounced hysteresis is observed in Fig. [21 for 5-values 
between 5.5 and 8.2. 

Fig. also contains the results of molecular dynam- 
ics simulations in a setup comparable to the experi- 
ment: K = 2 rectangular compartments with ground 
area 19.4 x 25.0 mm 2 and filled with TV = 400 particles 
of radius r = 1.18 mm and normal restitution coefficient 
v = 0.9 9] . While the slit at height h was as in the exper- 
iment, the hole at the bottom was chosen slightly larger 
(4.2 x 4.2 mm 2 ) to compensate for the vanishing thick- 
ness of the apertures and the neglected dissipation from 
particle- wall collisions in the simulations. Experiments 
and simulations are seen to agree well. 

These findings can be explained in the spirit of a pre- 
viously established flux model [7J with the help of two 
flux functions F and G describing the number of parti- 
cles which escape out of a compartment per time unit 
through the slit and the hole, respectively: 
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Here, B is defined in (JJ, n k is the fraction of parti- 
cles in compartment k £ {1,...,K} (hence n k £ [0,1], 
Sf^nfc = 1), and A essentially accounts for the size of 
the slit, [A] — s" 1 . The flux F(n k ) grows from zero for 




FIG. 2: Bifurcation diagram for the granular fountain. De- 
picted are the steady-state fractions of particles nh, k £ {1, 2}, 
in the two compartments from Fig. for different values of 
the driving frequency / in units of B oc /~ 2 from Eq. Q. 
Experimental results are indicated by the black error bars; 
the other symbols represent molecular dynamics simulations 
using uniform (blue asterisks) and clustered (red stars) initial 
conditions. The grey lines represent the theoretical predic- 
tion from the flux model for A = 0.018: solid lines 
correspond to stable configurations; dashed lines to unstable 



rih = to a maximum at n k = {BK 2 )^ 1 / 2 , and then de- 
creases asymptotically towards zero as a result of the in- 
elastic collisions between the particles. The flux through 
the hole, G(n k ), is obtained by taking the limit h — ► 
of F(n k ). The ratio A of G(n k ) and F(n k ) in the low 
density limit is determined by the ratio of the surface 
areas of the hole and the slit, and also accounts for the 
anisotropy of the velocity distribution just above the bot- 
tom. Obviously, A should be chosen considerably smaller 
than unity, otherwise the flux through the hole will com- 
pletely overpower that through the slit, and hence the 
clustering effect. For our experimental setup one finds as 
a rough estimate A « 0.02 under the assumption that the 
finite particle radius effectively reduces the dimensions 
of slit and hole by 1.8 mm and neglecting any velocity 
anisotropies. 

In the fountain geometry (K — 2) from Fig. the 
time-evolution of the relative populations n% = ni(t) and 
ri2 = 1 — ni follows by adding up all the fluxes through 
slit and hole: 



F(i - m) + G(i - m) - F(m) - G{ ni ) (3) 



The uniform distribution m = n-i = 1/2 is a stationary 
solution for all B, but is linearly stable only if the Jaco- 
bian-2F 1 '(l/2)-2G'(l/2) = -2A[(1-B) exp(-5)+A] is 
negative. This is the case when B is either small or large. 



If A < 



there exists a i?-interval for which the uniform 



distribution is unstable, and one can show that now the 
only stable solutions of (3J are a symmetric pair of clus- 
tered, fountain states. When fitting with A = 0.018, the 
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comparison of the predictions of the flux model with the 
experiment in Fig. [21 is fairly good; especially the bifur- 
cation and the hysteresis are well reproduced. 

Convection rolls somewhat similar to our present gran- 
ular fountain but governed by different physical mecha- 
nisms are well known for ordinary liquids and gases out 
of equilibrium and also for granular systems at very high 
densities. For a dilute granular gas, as we consider it 
here, convection in a single container has been reported 
recently |l0| , but neither a compartmentalized setup nor 
a theoretical modeling have been addressed before. 

We now proceed towards the "granular ratchet" setup 
in Fig. ^p: Starting with any number of copies of the 
granular fountain placed next to each other and closed 
cyclically, we then move each hole to its adjacent wall. 
Within the flux model © - © this amounts to replacing 
© by 

h k = G(n k+ i) + F(nfc_i) - F(n k ) - G(n k ), even k, 
h k =F(n k+ i) + G(n k -x)-F(n k )-G(n k ), odd k (4) 

with an even number of compartments K and periodic 
boundary conditions riK+ k = n k . Due to J2J) it follows 
that any periodically continued steady state solution of 
the original fountain geometry translates into a steady 
state solution of the extended granular ratchet geometry; 
the particle currents through the slits remain unchanged, 
and those through the holes are simply inverted. Less 
obvious is the stability of these solutions and the possible 
coexistence of further stable solutions: 

In the simplest case of a granular ratchet consisting of 
just K = 2 periodically closed compartments, one ob- 
tains, as expected, a bifurcation diagram (not shown) 
which practically coincides with that in Fig. [21 and sim- 
ilarly for the particle currents. Turning to K = 4 com- 
partments as depicted in Fig. we have performed 
molecular dynamics simulations with the same dimen- 
sions of compartment, hole, and slit and the same aver- 
age particle number N/K = 200 per compartment as in 
the K = 2 fountain case. The resulting bifurcation dia- 
gram in Fig. |3 displays as its most prominent new feature 
two different types of coexisting stable states with spon- 
taneously broken symmetry for not too large £>-values: 
(i) Ratchet states with alternating dense and dilute com- 
partments and a resulting finite directed particle flux, 
as quantitatively exemplified in Fig. 0] As the system 
evolves towards its stable steady state, both fluxes in 
Fig. 0] converge to the same constant mean. Its finite 
value is the most immediate signature of a spontaneous 
ratchet effect, (ii) Fluxless clustered states with alter- 
nating pairs of dense and dilute compartments and no 
resulting net particle flux, see Fig. 0] Our experimental 
data (not shown) on the K — 4 ratchet are in agreement 
with these findings. 

The intricate bifurcation diagrams for even larger K 
will be presented elsewhere. The main (and quite plau- 
sible) finding is that current carrying ratchet states with 
alternating dense and dilute compartments are always 
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FIG. 3: Bifurcation diagram of the K — 4 granular ratchet. 
Symbols: results from molecular dynamics simulations. Grey 
lines: prediction from the flux model (1), (2), (4) for A = 
0.018. Solid lines correspond to stable configurations; dashed 
lines to unstable ones. For further details regarding param- 
eters and methods see main text. For small B the uniform 
distribution, ni = 712 = nz — n& = 1/4, is the unique steady 
state. At B — 1 a pair of symmetry broken stable states takes 
over, namely 712 = ns ^ ru = m, indicated by blue asterisks. 
They exhibit no net particle flow and are called fluxless clus- 
tered states. For slightly larger B another pair of stable states 
emerges, both of the form m = 123 7^ 712 = n& (see Fig. Qd) , 
indicated by red stars. They do carry a net flow (cf. Fig. 
and are called ratchet states. For B > 6.8 only the fluxless 
stable states survive. 



stable solutions within an entire interval of i?-values. 
Moreover, one finds an increasing number of coexisting 
states with spontaneously broken symmetry. These are 
long-lived hybrid states, containing both ratchet-like re- 
gions (with alternating dense and dilute compartments) 
and fluxless clustered regions. 

Which specific steady state solution the system will 
settle in strongly depends on the initial condition, and at 
phase boundaries also on random fluctuations. An inter- 
esting example of the latter type is a uniform initial dis- 
tribution. The corresponding probability px of ending up 
in a ratchet state is depicted in Fig.El For an intuitive ex- 
planation we view the system as K/2 identical "dimers", 
each consisting of two compartments with a common slit 
at height h, and connected with each other by holes. 
Since the holes which connect neighboring dimers are 
much smaller than the slits, inside each dimer a quasi- 
steady distribution is approached much faster than the 
full equilibration between all dimers. Hence, an initially 
uniform distribution evolves inside each dimer towards 
either of two clustered states with probability q — 1/2, 
independently of the other dimers. This implies that a 
granular ratchet system with K compartments will end 
in one of the two ratchet states with alternating dense 
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FIG. 4: Running time average (over 5 seconds) of the net 
particle flux 4>(i) for the same molecular dynamics as in Fig. EH 
with B — 3.0, measured through the rightmost slit (a), and 
through the rightmost hole (b) in Fig. The upper (blue) 
curves in each plot correspond to the clustered state with 
densely populated leftmost and rightmost compartments in 
Fig. Et>. The lower (red) curves belong to the ratchet state 
with densely populated compartments 1 and 3 in Fig. 
For both states the straight line corresponds to the long-time 
asymptote of the particle current. 
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FIG. 5: Fraction pa of ratchet solutions versus number 
K (> 4) of compartments in a granular ratchet system with a 
constant initial particle number N/ K per compartment. Sym- 
bols represent numerical solutions of the (deterministic) flux 
model (1), (2), (4) for B — 3 and A = 0.018, averaged over ini- 
tial distributions with a superimposed randomization, mim- 
icking fluctuations due to the finite particle number within 
each compartment. The straight line is a linear best fit. 



and dilute compartments with probability 
PK = 2q*'\ 



(5) 



This prediction is in excellent agreement with the linear 
best fit p K = 2.02 (0.394) K / 2 in Fig. The reduced q- 
value can be attributed to the fact that the dimers are 
actually not entirely independent due the finite hole size. 

In practice, a ratchet state spontaneously arises only 
with a small probability for large K, as far as uniform ini- 
tial conditions are concerned. This situation is compara- 
ble to other systems with spontaneous symmetry break- 
ing, e.g. a ferromagnet [ll|: Just like a piece of iron 
picked from the shelf normally appears to be "unmagne- 
tized" and an external magnetic field is needed to pro- 
duce a magnetization on a global scale, a ratchet state 
only emerges with appreciable probability for large K if 



there is some initial bias towards this state. Without 
breaking the symmetry of the setup, this could be real- 
ized (in a reproducible direction, if so desired) by apply- 
ing a periodically modulated initial particle distribution 
or a small external force in the horizontal direction in 
Fig. acting during a certain preparatory time span. 

In conclusion, we have exploited the clustering effect 
to create spontaneous directed transport in symmetri- 
cally compartmentalized granular gases. In the granular 
fountain it takes the shape of a convection roll, and in 
the granular ratchet it appears as a current perpendic- 
ular to the vertical driving. In both cases, the directed 
transport arises as a collective effect of the stochastically 
colliding particles. 
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